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ABSTRACT 

We present observations and analysis of PSl-lObzj, a superluminous supernova (SLSN) discovered 
in the Pan-STARRS Medium Deep Survey at a redshift z = 0.650. Spectroscopically, PSl-lObzj is 
similar to the hydrogen-poor SLSNe 2005ap and SCP 06F6, though with a steeper rise and lower 
peak luminosity (Afb Q j ~ —21.4 mag) than previous events. We construct a bolometric light curve, 
and show that while PSl-lObzj's energetics were less extreme than previous events, its luminosity still 
cannot be explained by radioactive nickel decay alone. We explore both a magnetar spin-down and 
circumstellar interaction scenario and find that either can fit the data. PSl-lObzj is located in the 
Extended Chandra Deep Field South (ECDF-S) and the host galaxy is imaged in a number of surveys, 
including with the Hubble Space Telescope. The host is a compact dwarf galaxy (Mb w — 18 mag, 
diameter < 800 pc), with a low stellar mass (M* ~ 2.4 x 10 7 M Q ), young stellar population (t* « 5 
Myr), and a star formation rate of ~ 2 — 3 M Q yr _1 . The specific star formation rate is the highest 
seen in a SLSN host so far (~ 100 Gyr" 1 ). We detect the [O III] A4363 line, and find a low metallicity: 
12+(0/H) = 7.8 ± 0.2 (~ O.IZq). Together, this indicates that at least some of the progenitors of 
SLSNe come from young, low-metallicity populations. 
Subject headings: supernovae: general, supernovae: individual (PSl-lObzj) 



1. INTRODUCTION 

The discovery of "superluminous" supernovae 
(SLSNe), with peak luminosities 30 — 100 times 
brighter than normal supernovae and radiated energies 
> 10 51 ergs, is one of the most unexpected results 
from blank-field time-domain surveys like Pan-STARRS 
(PS1), Palomar Transient Factory (PTF) and the 
Catalina Real-Time Transient Survey (CRTS). Several 
distinct subclasses has been identified, indicating dif- 
ferent mechanisms to power the extreme luminosities. 
Some SLSNe can be classified as Type Iln, likely pow- 
ered by inter action with a dense, H-rich circumstellar 
medium (e.g. lOfek et all [20071: [Smith et al.l [20071 [2011 
iRest et all [201 It iMoriva et al.l 1201 21 ). The superlumi- 
nous SN2007bi was proposed to be a pair-instability 
explosion and so ult imately powered by radioactivity 
(|Gal-Yam et al.ll2009ft . though this claim is controversial 
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dDessart et all 12012ft. The rec ently discovered SLSN 
PSl-lOafx (jChornock et al.ll2013ft does not resemble any 
previous SLSNe and may define another class of objects. 

A third subclass of hydrogen -poor SLSNe s imilar 
to the transien ts SN2005ap (lOuimbv et al.l 12007ft 
and SCP 06F6 ([Barbarv et all [2009ft have~aIso been 
identified, characterized by blue spectra with a few 
broad featu res not matching a ny s t andard super- 
nova class (IQuimbv et al.l 120071 120111: IBarbarv et all 



20091: iPastorello et all 120101: iChomiuk et al.l 120111 
Leloudas et al.l 120121: iBerger et al.l 12012ft . While the 
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associated energetics, ejecta masses and host environ- 
ments point towards the explosion of a young, massive 
star, the ultimate energy source remains unknown 
for these objects. Like the H-rich SLSNe, models 
based on circumst el lar in t eraction has been pr o posed 
(IChevalier fc Irwinl I20T1 IGinzburg fc Balberd 120121: 
IMoriva fc Maedal 12012ft but the lack of hydro gen seen 
in the spectrum requires such interaction to be dom- 
inated by intermediate-mass elements. Alternatively, 
the luminosity could be explained by energy injection 
from a central engi ne, such as the spin-down of a 
newborn magnetar (iWooslevi 120101 : iKasen fc Bildstenl 
[20101 iDessart et all 12012ft . This class has also been 
linked to Type Ic SNe through the late-time spec- 
troscopic evolution of a few objects (jPastorello et al.l 
120101: iQufmbv et al.H2011ft . but the relationship between 
the classes remains unclear. Exploring the diversity 
of SLSNe and mapping the distribution of explosion 
properties will be important in further shedding light on 
the possible energy sources. 

Another clue to the progenitor systems could come 
from studying the host environments. Of the 10 
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2005ap-like H-poor SLSNe published p rior to this paper 
only five have detect ed host galaxie s (iNeill et alj I2011t 



Leloudas et "all l2012t iBerger et all I2012t iChen et al l 



2013( 1 Thehost galaxy of SN2010gx is the only one 



that has been studied in detail so far, and is a dwarf 
galaxy with a low metallicity [Z — 0.06Z Q ), leading to 
speculation of w hether metallicity plays a role in the pro- 
genitor channel (|Chen et al.ll2013f ). Increasing and char- 
acterizing the sample of SLSN host galaxies is essential 
for testing this hypothesis, and constraining the possible 
progenitors to these extreme explosions. 

Here, we present the discovery and analysis of PST 
lObzj, a hydrogen-poor SLSN at z = 0.650 from the Pan- 
STARRS Medium-Deep Survey (PS1/MDS). We present 
a comprehensive study of the SN and its host envi- 
ronment. The discovery and observations of PSl-lObzj 
are described in We analyze the properties of the 
supernova, including temperature evolution, bolometric 
light curve, possible models, and spectral modeling, in 
fj3j Since PSl-lObzj is located in the Extended Chandra 
Deep Field South (ECDF-S), its host galaxy i s detected 
in the GEMS, GaBoDs and MUSYC surveys (IRix et al.1 
12004 iTaylor et al.ll2009t iCardamone et al.ll2010fl . as well 
as in the PS1 pre-explosion images. The host galaxy 
properties, including metallicity, star formation rate, 
stellar mass and population age, are analyzed in 52J We 
place this SN in a broader context, comparing it to previ- 
ous events, and summarize our results in $5] and $6l All 
calculations in this paper assume a ACDM cosmology 
with H n = 70 km s" 1 Mpc" 1 , n M = 0.27 and fi A = 0.73 
(jKomatsu et al.ll2011[) . 

2. OBSERVATIONS 

2.1. PS1 Survey Summary 

The PS1 telescope on Haleakala is a high-etendue wide- 
field survey instrument with a 1.8-m diameter primary 
mirror and a 3.3° diameter field of view imaged by an ar- 
ray of sixt y 4800 x 4800 pixel detectors, with a pixe l scale 
of 0.258" (Kaiser et al.ll2010HTonrv fc Onakal l2009) . The 
observations are obtained through five broad-band filters 
(ffpiTpiZpi-zpiypi), with some differences relative to the 
Sloan Digital Sky Survey (SDSS); the g-pi filter extends 
200 A redward of gsDSS to achieve greater sensitivity and 
lower systematics for photometric rcdshifts, and the zpi 
filter terminates at 9300 A, u nlike zsnss which is de- 
fined by the detector response (jTonrv et al. 2012). PS1 
photometry is in the "natural" PS1 magnitude system, 
m = — 2.51og(i 7 ly) + m', with a single zero-point adjust- 
ment (m') in each band to conform to the AB magnitude 
scale, determined with PS1 ob servations of HST C alspec 
spectrophotometric standards (jBohlin et al.ll2001l ) . Mag- 
nitudes are interpreted as being at the top of the atmo- 
sphere, with 1.2 airmasses of atmospheric attenuation 
inclu ded in the system response function (jTonrv et al.l 
I20H . 

The PS1 Medium-Deep Survey (MDS) consists of 10 
fields (each with a single PS1 imager footprint) observed 
in gpirpiipizpiwith a typical cadence of 3 d in each fil- 
ter, to a 5cr depth of ~ 23.3 mag; j/piis observed near 
full moon with a typical depth of ~ 21.7 mag. The stan- 
dard reduction, astrometric solution, and stacking of the 
nightl y images is done by the Pan-STA RRSl IPP sys- 



tern (lMagnier| |2006; Magnic r" et al.l 120081 ) on a computer 



cluster at the Maui High Performance Computer Cen- 
ter. The nightly Medium Deep stacks are transferred 
to the Harvard FAS Research Computing cluster, where 
they are processed through a frame subtraction analy- 
sis using the photpipe pipeline devel oped for the Su- 
perMACHO and ESSENCE surveys J Rest et al.l 120051 
iGarg et al.l 120071: iMiknaitis efall 120071 ). which was fur- 
ther improved in order to increase the accuracy (Rest 
et al, in preparation, Scolnic et al., in preparation). The 
discovery and data presented here are from the photpipe 
analysis. 

2.2. Photometry 

PSl-lObzj was discovered in PS1 MD02 data on the rise 
on UT 2010 Dec 16, at coordinates RA=03 h 31 m 39.826 s , 
Dec=-27°47'42.17" (J2000). Spectroscopic follow-up 
confirmed it to be at redshift z — 0.650 from host galaxy 
emission lines, placing the peak observed absolute mag- 
nitude at < —21 mag, thus classifying it as "superlu- 
minous". The transient was detected in gpiTpiipi-zpi 
until PS1 stopped observing the field in early February 
2011. All photometry is listed in Table [1] and is cor- 
rected fo^OTe^rouna^rtiMtjpjowith E(B — V) = 0.008 
mag (jSchlaflv fc Finkbeinerl 120111 ). When PS1 resumed 
observing this field in late September 2011, PSl-lObzj 
had faded below the detection limit of ~ 23.5 mag. 

In addition to the PS1 photometry, griz images were 
obtained along with spectroscopic observations with the 
Low Dispersion Survey Spectrograph (LDSS3) on the 
6.5-m Magellan-Clay telescope, the Ina mori-Magellan 
Area ! Camera and Spectrograph (IMACS; Dressier et al. 
2006) on the 6.5-m Magellan-Baade telesco pe, and the 
Gem ini Multi-Object Spectrograph fGMOS; iHook et al.1 
2004) on the 8-m Gemini-South telescope, allowing 
us to extend the light curve until early April 2011. 
These images were reduced using standard routines in 
IRAFr^l. and transient flux was determined by subtract- 
ing the PS1 pre-explos ion template images using ISIS 
(jAlard fc Luptonl [1998) to correct for galaxy contami- 
nation, and measuring the flux in the difference image 
using aperture photometry. For the LDSS3 and IMACS 
images, calibrations were obtained either from observa- 
tions of standard fields on the same night, or from the 
PS1 catalogs of stars in the field of P Sl-lObzj corrected t o 
the SDSS system by the relations in lTonrv et al.l (|2012h . 
In the case of Gemini, archival zeropoints were used for 
calibration, after verifying that they produce consistent 
results with the PS1 catalog. 

In general the slight difference between the PS1 fil- 
ter set and griz would not introduce any significant er- 
rors. At the particular redshift of PSl-lObzj, however, 
the [O III] A5007 galaxy emission line is located at the 
edge between the i and z bands, contributing primarily 
to the zpi-filtcr in the PS1 photometric system, but to 
the i-band filter in the SDSS system used at Magellan 
and Gemini. This line contributes a substantial frac- 
tion of the galaxy flux (see SJ). Therefore, non-PSl i 
and z fluxes were either determined by subtracting the 
galaxy templates taken at Gemini and Magellan after 

10 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Observed light curve of PSl-lObzj. Time is shown in 
observer frame, relative to 2011 Jan 1 (MJD 55562.5). Circles show 
PS1 photometry, while the squares are photometry obtained with 
Magellan or Gemini. The arrows and triangles similarly show 3<r 
upper limits from PS1 and non-PSl photometry, respectively. We 
note the rapid rise time and the faster fall-off in the bluer bands, 
which indicates temperature evolution. 



the supernova had faded (§ 12. 4p . or corrected according 
to numerical subtraction. 

Figure [T] shows the observed light curves. Since PS1 
was observing this field prior to the detection, we are 
able to constrain the rise time, particularly in ipj-band, 
where PSl-lObzj brightened by > 1.2 mag in 12 days in 
the observed frame, corresponding to just 7 days in the 
rest frame. We also note that the later peak times in 
the redder bands indicate temperature evolution. Since 
the best-fit peak is different in different bands, we fit a 
low-order polynomial to our constructed bolometric light 
curve (§ 13. 2[) to determine the time of maximum light as 
UT 2011 Jan 02.65 (MJD 55563.65) ±2 d. All phases 
listed are in rest-frame days with respect to this zero- 
point. 

2.3. Spectroscopy 

We obtained four epochs of spectroscopy of PSl-lObzj. 
Details are given in Table [3J Our initial spectra were 
taken on 2011 Jan 18.2 using the Low Dispersion Sur- 
vey Spectrograph-3 (LDSS3) on the 6.5-m Magellan Clay 
telescope. Subsequent observations were obtained with 
the Gemini Multi-Object Spe ctrograph (GMOS ) on the 
8-m Gemini-South telescope (|Hook et al.l 120041 ). Con- 
tinuum and arc lamp exposures were obtained imme- 



diately after each object observation to provide a flat 
field and wavelength calibration. Basic two-dimensional 
image processing tasks were accomplished using stan- 
dard tasks in IRAF. Observations of spectrophotometric 
standard stars were obtained on the same night as the 
LDSS3 data, while archival observations were used for 
the GMOS spectra. Our own IDL routines were used to 
apply a flux calibration and correct for telluric absorp- 
tion bands. 

The LDSS3 observations covered the range 3540 — 
9450 A in a single setup using the VPH-all grating and 
a 0.75" slit oriented at the parallactic angle. Although 
no order-blocking filter was used for the object obser- 
vations, this setup exhibits very little second-order light 
contamination, which we confirmed from observations of 
standard stars taken both with and without a filter, so 
we believe the spectral shape to be reliable. 

The Jan 25 and 28 GMOS observations were taken 
with complementary blue and red setups, which we will 
sometimes present as a combined single spectrum. The 
Jan 25 blue spectra were taken with the slit oriented at 
a position angle of 175°, about 68° a way from the par- 
allactic angle, so differential light loss ([Filippcnko 119821) 
makes the blue continuum slope on that date unreliable. 
The other GMOS spectra were acquired in red setups ei- 
ther at low airmass (Jan 28) or near the parallactic angle 
(April 2 and 3), so their spectral slopes are reliable. 

Our last GMOS observations on Apr 2 
and 3 were obtained in no d- an d- shuffle mode 
(jGlazebrook fc Bland-Hawthorn! 1200 ID . An error 
resulted in the object being nodded off the slit for half 
of the Apr 2 observations. The exposure time quoted in 
Table [5] reflects only the on-slit time. The Apr 2 and 3 
data were combined into a single spectrum. 

The spectra from Jan 18, 25 and 28 are shown in Fig- 
ure [2] The Apr 2 spectrum is dominated by host galaxy 
light, and is shown in SjU All of our spectra show a 
number of narrow emission lines originating in the host 
galaxy, allowing us to determine a consistent redshift of 
z = 0.650 for PSl-lObzj. 

2.4. Host Galaxy Photometry 

The host galaxy of PSl-lObzj is detected in the 
PS1 pre-explosion stacked images in <?pirpiipi2pi. 
In addition, we obtained deep host images with 
Gemini-S/GMOS, Magellan-Clay/LDSS3 and Magellan- 
Baade/IMACS in griz after the supernova had faded. 
Deep infrared imaging in J and K with the FourStar In- 
frared Camera on Magellan-Baade (|Persson et al.ll2008f) 
yielded only upper limits. Table [3] lists all galaxy pho- 
tometry. 

PSl-lObzj is located in the Extended Chandra Deep 
Field South (ECDF-S) and so photometry from a num- 
ber of other surveys is also available. From the GEMS 
survey, there is Hubble Space Telescope (HST) imaging 
with the Adv anced Camera fo r Surveys (ACS) in F606W 
and F850LP (|Rix et all [20041 ). and we retrieved the re- 
duced images from the Mikulski Archive fo r Space Tele- 
scope s (MAST). From the GaBoDS survey (jTavlor et all 
2009]), there are detections in U^U BV RIz' (and non- 
detec tions in JHK). In ad dition, the MUSYC sur- 
vey (|Cardamone et al.l I2010D provides imaging in 18 
narrow-band filters. This field is also covered by the 
Spitzer Infrared Array Camera (IRAC) SIMPLE survey 
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(|Damen et al.ll20lTI ). but the galaxy is not detected in 
their catalog. Archival photometry of the host from the 
catalogs of these surveys is also included in Table |3l 

The various photometric measurements agree well in 
the bluer filters, but in the different i- and z-bands there 
is considerable discrepancy (e.g. the ipimeasurement is 
~ 0.6 mag fainter than the GMOS i-band, while the 
zpi measurement is ~ 0.9 mag brighter than the corre- 
sponding z filter). This is explained by the redshifted 
[O III] A5007 emission line, located near the edge be- 
tween i and z. The flux we measure in this line from 
the spectra (Sj4| is consistent with the differences in pho- 
tometry. The effect of this line is also clearly seen in the 
narrow-band photometry in the IA827 filter. 

3. SUPERNOVA PROPERTIES 
3.1. Spectroscopic and Light Curve Comparisons 

Given the redshift of z = 0.650, we find that PST 
lObzj reached a peak absolute magnitude of —21.17 ± 
0.15 mag in gpi. This is luminous enough to be 
classified as "supe rluminous " acco rding to the defini- 
tion sugg ested in iGal-Yaml (|2012f ). Figure H shows 
our spectra of PS l-lObzj, compared to hydrogen-poor 
SLSNe PSl-lOky (IChomiuk et alJlMLll) an d SN2010gx 
(jPastorello et alll201ft iQuimbv et all 120111 ). Our first 
spectrum shows a blue continuum with broad UV fea- 
tures that a re characteristic of t he class of hydrogen - 
poor SLSNe (jQuimbv et al.ll201ll : IChomiuk et al.ll20lTI ). 
These features are also visible in the blue GMOS spec- 
trum taken 7 rest-frame days later. In the red GMOS 
spectrum on day 16, a number of broad, low amplitude 
features have also developed, similar to the features seen 
in SN2010gx at late time. The combination of its lumi- 
nosity and spectral features unambiguously establishes 
PSl-lObzj as another member of the class of 2005ap-like, 
hydrogen-poor SLSNe. 

Figure [3] shows the light curve of PSl-lObzj in absolute 
magnitude versus rest frame phase compared to a few 
other hydrogen- p oor S LSNe: PSl-lOky and PSl-lOawh 
(IChomiuk et all 120111) SN2010gx dPastorello et al 
2010t IQuimbv et al.ll201lD and PTF09cnd~ i|Quimbv et al 
201 If) . We do not carry out detailed k-corrections due 
to the uncertainties in the spectral energy distributions 
(SEDs), but have picked bands at similar effective wave- 
lengths as indicated on the plots to facilitate compar- 
isons. With a fast rise time and slower decline, PSl-lObzj 
does not show the clearly symmetric ligh t curve behav- 
ior se en in previous hydrogen-poor SLSNe (jQuimbv et all 
120 111) , though we note that the rise time is less well con- 
strained in the bluer bands due to shallower limits prior 
to detection. In general, the light curve of PSl-lObzj ex- 
hibits similar timescales to SN2010gx and PSl-lOky, but 
has a flatter peak and is fainter overall. 

3.2. Temperature Evolution and Bolometric Light Curve 

We determine blackbody temperatures by fitting 
Planck functions to the broadband photometry, using a 
^-minimization procedure. For the PS1 photometry, 
where different bands are observed on consecutive rather 
than the same night, we first interpolate the photome- 
try to a common time by fitting a low-order polynomial 
to the nearby light curve points. The spectral energy 
distribution (SED) fits are shown in Figure |H with the 
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Figure 2. Spectra of PSl-lObzj (black), c ompared to spectra o f 
other hydrogen-poor SLSNe PSl-lOky (blue ; IChomiuk eTaT][20Tll) 
and SN2010gx (red; IPastorello et al .1 [20101') . The blue continuum 
and broad UV features are common to hydrogen-poor SLSNe. 
By the Day 18 spectrum, PSl-lObzj had also developed a num- 
ber of features in the optical , sim ilar to those seen in SN 2010gx 
(Pastorcllo et al. 2010). Sec i|3.3l and for modeling and identifica- 
tion of the features. 
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Figure 3. Absolute magnitude light curve of PSl-lObzj (solid 
black circles and black arrows) compared to other hydrogen- 
poor SLSNe at similar rest-frame wavelengths, showing gpi(~ 
3000 A rest-frame, top), rpi(~ 3800 A rest-frame, middle), and 
ipi(~ 4600 A rest-frame, bottom). Blue diamonds and trian- 
gles show PSl-lOawh and PSl-lOky rpi, ipiand zpi respectively 
(Chomi uk ct al. 2011), red squares show SN 2010gx in it, g and r 
(Pastorcllo ct al. 2010; Quimby et al. 2011), and yellow open cir- 
cles show PTF09cnd at ~ 3600 A dOuimbv et alj|201ll ). 
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model temperatures and radii indicated. These numbers 
should be interpreted with some caution, as especially 
at later times the spectrum clearly deviates from that 
of a blackbody. In addition, by the time of the first 
spectrum the broad UV absorption features is clearly af- 
fecting the g-band flux, so that the temperature inferred 
from photometry is lower than that found by fitting to 
the spectrum. 

The resulting blackbody temperatures and radii from 
all the fits to the photometry is shown in the top two pan- 
els of Figure [5j Prior to peak, we can only place a lower 
limit on the temperature of ~ 15,000 K, since the peak 
of the blackbody curve is bluewards of the observed pho- 
tometry and we are essentially fitting the Rayleigh- Jeans 
tail. After the peak, we find a clear decline in tempera- 
ture and increase in radius. The best-fit straight line to 
the estimated blackbody radius (Figure[5]) corresponds to 
an expansion velocity of 11, 000 ± 2000 km s _1 , in good 
agreement with velocities derived from spectroscopic fea- 
tures (§E2]). We note that dividing the estimated radius 
at peak by the velocity gives a timescale of ~ 23 days, 
consistent with the observed photometric rise. 

To construct a bolometric light curve, we first sum the 
observed flux by trapezoidal integration, interpolating to 
the edges of the observed bands. Since this only takes 
into account the flux in the observed wavelength range, 
it should be considered a strict lower limit of the total 
radiated power. Integrating this luminosity over the time 
period we observed the SN indicates a lower limit of the 
radiated energy E, ad > (2.4 ± 0.5) x 10 50 ergs. The 
resulting light curve is shown as open red circles in the 
bottom panel of Figure [5] 

To improve this estimate, following iChomiuk et al.1 
(|2011[) we also include a blackbody tail redwards of the 
observed bands, using the temperatures determined by 
our blackbody fits. The resulting pseudo-bolometric light 
curve is shown as the black, filled circles the bottom panel 
of Figure [5j We also include the early ipidetection, as- 
suming the same bolometric correction as the next light 
curve point. PSl-lObzj reached a peak bolometric mag- 
nitude Mbol = —21.4 ± 0.2 mag, and an estimated total 
radiated energy of i? ra d > (3.5 ± 0.6) x 10 50 ergs. As 
expected from the light curve, this is significantly less 
luminous than previous events - for example, PSl-lOawh 
and PSl-lOky reached peak bolometric magnitudes of 
—22. 2 mag and —22.5 mag respectively (Chomiuk ct al. 
1201 If) . While the spectroscopic features clearly identifiy 
PSl-lObzj as a member of the same class of objects, it is 
one of the least luminous hydrogen-poor SLSNe discov- 
ered to date. 

3.3. Line Identifications 

We used the supernova spectrum synthesis code SYNOW 
to obtain line identifications and estimates of the ex- 
pansion velocities, including manual and automated pro- 
cedures employing the recently updated versions of the 
software SYN++ in combination with SYNAPPSH The ba- 
sic assumptions of SYNOW include spherical symmetry, ve- 
locity proportional to radius, a sharp photosphere, line 
formation by resonant scattering treated in the Sobolev 
approximation, local thermodynamic equilibrium (LTE) 
for the level populations, no continuum absorption, pure 

1 Software was retrieved from https://c3.lbl.gov/es/ 
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Figure 4. Spectral Energy Distribution fits to the photometry. 
The best-fit blackbody temperatures and radii are indicated in the 
individual panels. The concurrent spectrum is shown with the 
photometry in the fourth panel; we adopt the temperature derived 
from the spectrum rather than from the photometry for this date. 

resonance scattering, and only thermal excitations. Fits 
are constrained by how we are able to best match absorp- 
tion minimum profiles, as well as the relat ive strengths 
of all the features (see Br anch et alj 2002 fo r more de- 
scription of fitting parameters and [Thomas et alJ 1 2 1 ll 
for software details). 

The fit to the 2011 Jan 13 spectrum (phase +7 d) is 
shown in the left panel of Figure [6] The photospheric 
velocity is set at 13,000 km s _1 , and the temperature to 
15,000 K. A maximum cut-off velocity of 40,000 km s _1 
was used for all ions, with the minimum velocity set to 
13,000 km s _1 . These parameters are comparable to the 
ions and associated veloci t ies identified in other ULSNc 
(e.g.. iQuimbv et all [20lH: IChomiuk et~aH IMtl . Two 
strong features observed around 2440 and 2650 A are 
fit reasonably well with Si III and Mg II, respectively. 
Introduction of Fe II improves the fit around the Si III 
line, as seen in the inset. Without Fe II, the red wing of 
the absorption could not be fit with Si III alone. 

Additional weaker features with less certain identifica- 
tions are also seen. A sharp cut-off around 3000 A is 
likely attributable in part to Si III, and we fit absorption 
features around 4230 and 4490 with O II. We include in 
the synthetic spectrum C II, which is cut off to the blue 
of Si III, but is seen in the other SLSNe and in the later 
spectrum of this object. 
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Figure 5. Top panel: Temperature evolution of PSl-lObzj, as 
determined from fitting a blackbody curve to the observed pho- 
tometry. The uncertainty at early times is largely due to the peak 
of the blackbody being bluewards of our bluest bands; see Fig- 
ure [4] Middle panel: Radius of PSl-lObzj, as measured from the 
same blackbody fit to photometry as the temperature. The best- 
fit straight line (dashed) corresponds to an expansion velocity of 
11, 000 ±2000 km s _1 . Bottom panel: Estimated bolometric light 
curve of PSl-lObzj. The open red circles show observed flux only, 
while the black filled circles include the observed flux plus a black- 
body tail in the red. 

The 2011 Jan 25 and 28 spectra were combined to fit 
a single phase +15 d spectrum, shown in the right panel 
of Figure [6] The photospheric velocity is set at 11,000 
km s _1 and the temperature to 11,500 K. The maximum 
cut-off velocity was once again set to 40,000 km s -1 for all 
ions, which were fitted with minimum velocities ranging 
between 11, 000 — 15, 000 km s _1 . We observe the same 
absorption features associated with Si III, Mg II, and 
C II. Again, including Fe II substantially improves the 
fit in this region; the contribution from Fe II only to the 
fit is shown in the inset. However, the lack of additional 
lines elsewhere in the spectrum at this temperature, as 
well as the noise at the bluest wavelengths, prevent de- 
termining the relative strengths accurately. Additional 
weaker features are fit with Ca II and Si II. The O II 
seen in the earlier spectrum no longer appears to be a 
conspicuous contributor to the spectrum. 

Two other SLSNe have shown significant spectral evo- 
lution post-peak : in both SN2010gx and PTF09cnd 
(jPastorello et al.l 120101 : iQuimbv et all 120111 ). the spec- 
tra evolved to look like normal Type Ic SNe at late 
times. The link to SNe Ibc is also suggested by 



the presence of similar, transient O II features in 
SN2008D, a well-observed SN lb (e.g.. iSoderberg et al.l 
l2lMlModiaz et alj|2009t iPastorello et al.ll2010l) . A basic 
question is whether the SLSNe are truly distinct objects 
from normal Type Ic SNe, or whether there is a smooth 
continuum between the two. With its comparatively low 
peak luminosity, PSl-lObzj is closer in luminosity to lu- 
minous Type Ic SNe like SN 20 1 Pay, which peaked at 
M r = -20.2 ([Sanders et al.ll2012[ ). than to the prototype 
SLSNe 2005ap and SCP 06F6. From this perspective, it 
is interesting to note that at least over the time we were 
following it, the spectral features in PSl-lObzj do not 
look like SN Ic features, including objects like SN 2010ay. 

3.4. Light Curve Model fits 

The optical luminosity of most canonical Type I SNe 
(i.e. type la, lb and Ic) is powered by the radioactive 
decay of 56 Ni, with the shape of the light curve primar- 
ily dictated by three parameters: the nickel mass (Mni) 
which sets the total luminosity, the total kinetic energy 
(Ek), and ejecta mass M e j, which set the characteris- 
tic time of photon diffusion r c oc M^E^ 1 ^ 4 and essen- 
tially determines the width of the light curve ([Arnettl 
|1982( K Measurements of the photospheric velocity (v p h) 
from the spectra constrain UE^JM^, so that all three 
parameters can be determined based on observable quan- 
tities. We fit our bolometric light curv e of P Sl-lObzj us- 
ing th e models of iValenti et al.l (|2008l ) and iDrout et al.l 
( 2011); see Figure [7] The light curve can be reason- 
ably fit with MNi — 6 — 8 Mq, with the best-fit model 
having Mm — 7.2 M and r c ~ 19 d. Using the pho- 
tospheric velocity derived from the spectrum near peak, 
«ph = 13,000 km s _1 , yields M ci - 5 - 11 M , with 
8.5 M Q for the best-fit model. Therefore, if PSl-lObzj 
were powered by radioactive decay, it would require a 
56 Ni mass of > 10 times what is o bserved in typ ical Type 
Ibc or Ic-BL SNe (0.2 - 0.5 M Q ; IDrout et alJlMTl . In 
addition, the ejecta would have to be 75 — 100% 56 Ni 
by mass, a fraction seen in no observed SNe, including 
proposed pair-instability SNe lik e SN2007bi where th e 
inferred Ni mass was several M (|Gal-Yam et al.ll2009f ). 
A composition of > 75% Ni would also result in a large 
amount of line-blanketing in the rest-frame UV, which 
is not seen. Radioactive decay, then, is unlikely to be 
the main contributor to the luminosity. This is consis- 
tent with what is found for other hydroge n-poor SLSNe 
(jChomiuk et al.ll27)TTI : IQuimbv et alll201ll ) . 

Since nickel decay is unlikely, other explanations have 
been proposed for the extreme luminosity of SLSNe. One 
possibility is energy injection by a cent ral engine, such as 
the sp in-down of a newborn magnetar (jKasen fe Bildstenl 

EMa |^Qsk3 [io!3L_We_fit. our bolomet ric light curve 

with the model ofK asen fe Bildstenl (I20101). fo llowing the 
procedure outlined in Chomiuk et al.1 (|201lD . Our as- 
sumptions include magnetic dipole spin-down, an opac- 
ity of k = 0.2 cm 2 g -1 , and a supernova energy of 
10 51 ergs; we vary the ejecta mass, the magnetar spin 
(p), and the magnetic field (B). With these assump- 
tions, Figure [8] shows a best-fit model, with M C j — 2 M , 
B = 3.5 x 10 14 G, and p = 4 ms; however we find that the 
light curve can be reasonably fit within the uncertainties 
with ejecta masses in the range M G j ^1 — 6 Mq. An 
ejecta mass lower than 1 M Q predicts a light curve that 
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Figure 6. SYNOW fits to the spectra of PSl-lObzj. Left: The 2011 Jan 13 (phase +7 d) spectrum. The strong UV features are well 
fit with Si III and Mg II; minor features near 4500A could be due to O II. The shape of the Si III feature is better fit when including 
Fe II, as illustrated in the inset, which shows the fit without including this ion. Right: Fit to combined 2011 Jan 25+28 (phase +15 d) 
spectrum. Minor features in the optical are fit with Ca II, Mg II and Si II. The fit includes a substantial contribution from Fe II (shown in 
the inset), but the overlap of these features with the Mg and Si features and lack of features elsewhere prevent us from determining the 
relative strengths accurately. 
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Figur e 7. Mod el of radioactive Nickel decay, following 
IDrout et al. (2011). The light curve can be reasonably fit with 
Ni masses in the range 6 — 8 Mq, but require the ejecta composi- 
tion to be 75 — 100% Ni to simultaniously fit the peak luminosity 
and the light curve width. 

is too narrow, while ejecta masses greater than ~ 6 M© 
require initial spins faster than the maximum (breakup) 
spin of ~ 1 ms to match the timescales. Within this 
range, however, parameters can be chosen to fit the light 
curve equally well within the uncertainties. We note that 
the inferred ejecta ma sses are similar to what is seen in 
normal Type Ibc SNe (jDrout et al.ll2011l ). 

The magnetar model predicts that the ejecta will be 
swept up into a dense shell, which then expands at a con- 
stant velocity. For the range of models that fit the light 
curve, those with higher spin periods (and so a lower total 
energy, as the the magnetar contribution scales as p~ 2 ) 
generally have lower inferred velocities, in better agree- 
ment with the velocities inferred by the spectra. The 
2 M© model shown in Figure [8] has the swept-up shell 
expanding at 11,000 km s _1 , in good agreement with the 
observed velocities. The predicted temperatures for this 
model also match the observed temperatures within the 
errors. A simple magnetar model thus provides a rea- 
sonable fit to the the observed properties of PSl-lObzj. 
One caveat is that our modeling of the spectra support 
a slightly declining, rather than constant photospheric 
velocity. The rapid evolution of the spectrum is also 
challenging to explain in the context of this model. 

A third proposed mechanism for powering SLSNe is 
interaction with opaque, circumstellar material. This 



Figure 8. Magnetar model fit to the bolometric light curve, with 
M cj = 2 M , B = 3.5 X 10 14 G, and p = 4 ms. While the 
observed light curve can be fit with a range of parameters within 
the uncertainties, the model shown also predicts a velocity and 
temperature evolution that agrees with the observed data. 

leads to efficient conversion of the kinetic energy to 
radiation, with the resulting lightcurve being due to 
shock breakout through this opaque wind. This class 
of models has been applied both to superluminous 
SNe Iln such as SN2006gy, and t o SN 2005ap- 
like objects (ISmith fe McCrayl 120071: ISmith et all 
pOlOt IChevalier fe irwinl 1201 U IBalberg fc Loebl 120111: 
IGinzburg fe Balberd I2012t IChatzopoulos et al.l 12012ft . 
The light curve of the hydrogen-poor and superluminous 
SN2006oz, in particular, showed a "dip" feature on the 
rise that has been interp reted as a signature of shock 
breakout (|Leloudas et al.ll2012t iMoriva fe Maedall2012D . 

We can use the observ ed properties of PSl-lOb zj and 
the analytical relations of lChevalier fe Irwinl (|2011D to es- 
timate the physical conditions required in the interaction 
scenario. Assuming a wind density profile p w = Dr~ 2 , 
as expected from a steady wind, so that D = M/4irv w = 
5 x 10 16 £>* in cgs units. The rise-time can be roughly 
equated to the diffusion time, td = 6.6k£>* d, where k is 
the opacity in units of 0.34 cm 2 g" 1 . We use k — 0.59, 
as expected for an ionized He-rich wind. Taking the rise- 
time of PSl-lObzj to be ~ 20 days, we find that £>* ~ 5.1, 
which gives a total required wind mass of ~ 3.5 Mq, 
using the radius at peak to be ~ 2.2 x 10 15 cm. Us- 
ing E xa d « 3.5 x 10 50 erg, we find that the associated 
supernova energy is 2.2 x 10 51 (Af cj /10M Q ) 1 / 2 erg , and 



corresponding diffusion radius Rd ~ 1 x 10 15 cm, as- 
suming an ejecta mass of 10 and using Equations 5 
and 3 in lChevalier fc Irwinl (|2011l ) respectiveiy. The pre- 
dicted velocity of the photosphere, using Equation 3 in 
IChomiuk et al.l ([201 ID . is 11,800 km s _1 , in good agree- 
ment with the observed velocities. Thus, this model can 
also reproduce the basic observed properties, but require 
a wind mas s of several M ffi of hydrogen -poor material. 

Recently, iGinzburg fc Balberd (I2012T) have show n that 
the simple treatment in Chevali er fc Irwinl fj201 lh is not 
appropriate in regimes where the wind radius is com- 
parable to the diffusion radius. Instead, they carried 
out hydrodynamical simulations of supernovae exploding 
into dense circumstellar material, successfully matching 
the light curves of SN 2005ap, SN 2006gy and SN 2010gx. 
Given the similarities between PSl-lObzj and SN2010gx, 
it seems plausible that its light curve could also be fit 
by a more sophisticated shock breakout model, though 
calculating such a mo del is outside the scope of t his pa- 
per. We note that the IGinzburg fc Balberd (|2012f ) model 
for SN2010gx requires an even larger total wind mass 
(M w ~ 16 M ff) than our estima t e for PSl-lObzj based 
on the simple IChevalier fc Irwinl ([20111 ) relations, so an 
extreme mass loss episode would likely still be required. 

A simple interaction model, then, can also explain 
the observed data, but requires a mass-loss rate of ~ 
3 Mq yr -1 in the last year before explosion, assuming 
a wind vel ocity of 1,000 km s -1 (as seen in Wolf-Rayet 
stars; e.g. Nugis fc Lamersl 120001 ) . In addition, the lack 
of hydrogen and helium seen in the spectra would require 
this circumstellar material to be primarily composed of 
intermediate-mass elements. One might also expect to 
see intermediate- width lines in the spectra if the primary 
energy source is interaction, but this has not been seen 
in any of the H-poor SLSNe, including PSl-lObzj. A de- 
tailed radiative transfer model is necessary to see if this 
scenario can reproduce the spectra as well as the light 
curves. 

4. HOST GALAXY PROPERTIES 

In addition to studying the SN explosion itself, addi- 
tional clues to the nature of the progenitors come from 
studying the host environments of the SLSNe. In the 
case of PSl-lObzj there is a wealth of data on the host 
galaxy, allowing for a detailed study. 

4.1. Luminosity and Size 

The absolute magnitude of the host is Mb = 
— 18.0 mag, corrected for cosmological expansion and 
foreground extinction. This is similar to what has 
been seen for other SLSN hosts, whic h seem to show a 
preference for low -luminosity galaxies (jNeill et al.l 1201 It : 
IChen et al.l 20131). In terms of the luminosity fun ction 



at z ~ 0.7 (jllbert et al.ll2005t IWfllmer et al.l [20061. this 
corresponds to a 0.05L* galaxy. 

The host galaxy of PSl-lObzj is unresolved in all our 
ground-based images (with seeing down to ~ 0.6"). Since 
the field was covered by the GEMS survey, we also 
have availa ble HST/ACS images in F606W and F850LP 
(|Rix et al.l 120041 ) , shown in Figure |9] Even in these im- 
ages, the host is not obviously resolved with a FWHM of 
~ 0.12" (whereas the mean FWHM of stars in the images 
is ~ 0.10"). At z = 0.650, this corresponds to an upper 
limit on the galaxy diameter of < 800 pc. We note that 




Figure 9. HST/ACS images of the field around PSl-lObzj, in 
filters F606W (left) and F8501p (right). The host, marked with 
the blue arrow, is remarkably compact. 



it is possible that what we see in the HST images is only 
one bright knot of star formation, and the galaxy itself 
could be more extended. Nevertheless, the combinaton 
of luminosity and size establishes the host as a compact 
dwarf galaxy. 

4.2. Stellar Mass and Population Age 

To determine the stellar mass (-M*) a nd population age 
(r*) of the host we fit the SED with the lMarastonl (|2005[ ) 
evolutionary stellar population synthesis models, using a 
Salpeter initial mass function and a red horizontal branch 
morphology. Since the model only accounts for the con- 
tinuum emission, and the flux in the host emission lines 
is substantial, we restrict our fit to the MUSYC narrow- 
band filters without significant emission lines. The fit to 
the SED is shown in Figure [KjJ We find that the host 
SED is well fit with a stellar population age of t* ps 5 
Myr, yielding a stellar mass of M* 2.4 x 10 7 M Q . This 
assumes Ay = 0, measured from the Balmer decrement 



Another estimate of the stellar population age comes 
from the H/3 equivalent width (EW). While we do not 
have a galaxy-only spectrum, the 2011 Apr 3 Gemini 
spectrum (Figure fTT|l does not show any broad super- 
nova features and is dominated by galaxy light (> 50%, 
estimated from pre-explosion galaxy photometry). It can 
therefore be used to determine a lower limit on the Hf3 
EW, which we find to be W r ~ 61 A. This value yields a 
young stellar population age of < 5 Myr for a metallicity 
Z = 0.2-0.4 Zm . using the fits inlLeyesaue et all (|2010al) 
to the models of lSchaerer fc Vaccal (|1998l ). This value is 
in excellent agreement with the stellar population age 
inferred from SED modeling. 

4.3. Metallicity 

While all of our spectra include contributions from 
both the galaxy and the SN, they clearly exhibit narrow 
emission lines originating in the host galaxy. Figure 1111 
shows the 2011 April 3 GMOS spectrum, which is dom- 
inated by galaxy light, with the strongest emission lines 
marked. We measure the line fluxes in all of our spectra 
by fitting Gaussian profiles (Table [4j . With the excep- 
tion of the [O III] A4363 line, which was only robustly 
detected in the 2011 Jan 28 GMOS spectrum, we use the 
weighted average of the three measurements for line di- 
agnostics. Absolute flux calibration is based on the 2011 
Jan 13 LDSS3 spectrum, by scaling synthetic photome- 
try from the spectrum to photometry obtained the same 
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night. The GMOS spectra were then scaled according to 
the flux in the [O III] doublet. 

None of our spectra cover Ha, which is located at 
1.083 /jm at this redshift. We therefore use the Balmer 
decrement as measured from Hy/H/3 to estimate red- 
dening, assuming a Case B recombination value of 0.469 
(|Osterbrocld H989). Since our measured value of 0.48 ± 
0.03 is consistent with no reddening, we conclude that 
the host galaxy extinction is minimal. 

We detect the auroral [O III] A4363 line in the 2011 
Jan 28 GMOS spectrum at 3.5c significance, shown in 
the inset of Figure [TTJ Assuming an electron density 
n e = 100 cm -3 , we calculate an electron temperature of 
T e (0++) = 16, 2001^70° K from the ratio of [O III] A4363 
to [O III] AA5007,4959, using the IRAF task temden. 
This result is not sensitive to the exact choice of den- 
sity since T e is inse nsitive to small changes in density 
(|Kewley et al.1 120071 ) : for example, we find identical re- 
sults when doubling the assumed electron density to 200 
cm" 3 Using the rela tion T e (0+) = 0.7 x T e {0 ++ ) + 0.3 
from lStasiriskal (1 1982]). we determine Q+ /H and ++ /H 
using the relations in lShi et al.l ((2006) . This gives an elec- 
tron temperature metallicity of 12+log(0/H) = 7.8±0.2. 
This tra nslates to Z = 0A3Z Q , using the solar abun- 
dance of lAsplund et al.l (|2009D . This low abundance is 
consistent with the inferred young stellar population age 
and low stellar mass. 

For comparison, we also estimate the oxygen abun- 
dance using the i?23 diagno s tic w ith the calibra- 
tion of IKobulnickv fc Kewlevl (|200l . We measure 
R 23 = ([O III] AA5007, 4959 + [O II]A3727)/H/3 = 
9.25 ± 0.36, and an ionization parameter of y = 
log ([O III]AA5007, 4959/J O II]A3727) = 0.95 ± 0.03. Us - 
ing the iterative scheme in IKobulnickv fc Kewlevl (|2004|) . 
and assuming the lower metallicity branch based on the 
presence of the [O III] A4363 line, this method gives a 
metallicity 12 + log(0/H) « 8.3. This is 0.5 dex higher 
than the result from the direct method, but we note 
that this discrepancy is not unusual; theoretical strong- 
line indicators are known to be offset from the direct 
method, with t he difference being larger at the lower- 
metallicity end (|Bresolin et al.1120 09 1 ). A similar discrep- 
ancy is se e n in the host galaxy of SN2010gx, where 
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Stoll et alJ (I2011D foun d 12 + log(0/H) = 8.36 using the 
Kobulnickv fc Kewlevl (|2004l ) calibration, w hile the di- 
rect method yields 12 + log(0/H) = 7.46 (|Chen et all 
l2013f h Still, strong-line metallicity indicators provide a 
useful basis for comparison, since direct metallicity indi- 
cators are otherwise mostly only available for low-redshift 
samples. 

4.4. Star Formation Rate 

We estimate the star formation rate (SFR) of the 
host galaxy from the [O II] A372 7 line flux, us i ng th e 
metallicity-dependent relation in iKewlev et al.l (2004). 
Using the metallicity and ionization parameter we deter- 
mined from the R23 method, we find SFR « 2 M yr" 1 . 
Alternatively, since the H7/H/3 ratio indicates no extinc- 
tion, we can use the H/3 flux to predict the Ha flux, 
assuming a ratio Ha/H/3 = 2.85 according to case B 
recombination. U sing SFR = 7.9 x 10~ 42 Lna(erg s~ ) 
(|Kennicuttill99l . we find SFR w 4.2 M yr , in rea- 
sonable agreement with the [O II] A3727 estimate. 
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Figure 10. Best-fit spectral energy distribution models to the 
host galaxy photometry, using Ay = mag. Only the narrow-band 
photometry is used in the fit (filled circles) , due to the strong emis- 
sion lines contaminating the broadband niters. Open circles show 
narrow-band filters which contain emission lines and are therefore 
not used for the fit. For comparison, the stars show PS1 broad- 
band photometry, both uncorrected (open) and corrected (filled) 
for the flux in the emission lines. Additional broad-band filters are 
largely redundant with the PS1 ones and are not shown. Our best- 
fit is a young (r» m 5 Myr) and low-mass (M* « 2.4 x 10 7 M Q ) 
stellar population; an acceptable fit also exist for a 3.5 Myr pop- 
ulation, and slightly worse fits for a 10 or 20 Myr population, as 
shown in the inset. 



A complementary method to calculating the SFR is to 
use the galaxy UV continuum flux. At this redshift, the 
UBg filters sample rest-fra me 2300 — 2900 A , allowing 
us to use the relation from iKennicutti (|1998l ): SFR = 
1.4 x 10~ 28 L„. This yields SFR w 2 - 3 M yr" 1 , also in 
good agreement with the estimates from emission lines. 

Combining the stellar mass with the star formation 
rate, we calculate a specific star formation rate (sSFR) of 
~ 100 Gyr -1 . This is significantly higher th an the ~ 2.6 
Gyr -1 measured in the host of SN 2010gx (|Chen et al.1 
I2013D . and also higher than the ~ 10 Gyr -1 in the host 
of PSl-lObam (|Berger et al.l 1201 21 ) . The basic picture of 
a low metallicity, low mass and highly starforming dwarf 
galaxy is similar to what has been seen for other SLSNe. 



5. DISCUSSION 



5.1. 



The Diversity of SLSNe 

As was shown bv lQuimbv et al.1 (|2011[) . the hydrogen- 
poor SLSNe form a spectroscopic class, though with a 
range of light curve properties. Figure [12] shows the dis- 
tribution of peak absolute magnitudes of all published 
2005ap-like hydrogen-poor SLSNe, corrected for cosmo- 
logical expansion by M — m — 5 log (dr,(z)/10pc) + 
2.51og(l + z). Due to the lack of SED information in 
several objects, we do not carry out full k-corrections, 
but note that where multiband photometry is avail- 
able the observed peak is at a rest-frame wavelength 
of - 3000 - 4000A. SN 2006oz is not included in this 
plot, as it was only observed on the rise and so the 
peak magnitude is not well constrained. Most of the 
hydrogen-poor SLSNe peak near M ~ — 22 mag, with a 
tail to higher luminosities. The apparent lack of lower- 
luminosity objects is likely due, at least in part, to the 
flux-limited surveys (and spectroscopic follow-up) so that 
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Figure 11. 2011 Apr 3 CMOS spectrum of PSl-lObzj, with strong 
host galaxy emission lines marked. While there is some supernova 
contribution to the flux in this spectrum, it is dominated by host 
galaxy light (> 50%). The inset shows the region around the au- 
roral [O III] A4363 line, from the 2011 Jan 28 spectrum, which has 
somewhat better S/N. 

there is a bias towards finding brighter objects. PSl- 
lObzj is both the lowest-luminosity and one of the lowest 
redshif t SLSNe found in P S1/MDS, for example. Re- 
cently, [Quimb^eFal- (2013) found that the distribution 
of hydrogen-poor SLSNe seems to be narrowly peaked, 
also when taking the effects of flux-limited selection into 
account. If so, an event like PSl-lObzj would be in- 
trinsically rarer than the —22 mag objects, at the low- 
luminosity tail of the distribution. 

The timescales seen in SLSNe also vary by a factor of 
several, with rise times varying from ~ 2 d in the case 
of PS l-lObzj, to > 50 d in PTF09cnd (|Quimbv et al.l 
2011). If the faster timescales are typical for the lower- 
luminosity end of the distribution, it may present an ad- 
ditional selection bias against the fainter objects, as the 
timescales are approaching those of normal SNe and so 
the objects stand out less amongst the more common 
normal SNe. 

Finding fainter hydrogen-poor SLSNe is particularly 
interesting because this class is linked to Type Ic SNe 
through the late time spectroscopic evolution of a few 
objects. A basic question is whether they are truly dis- 
tinct populations, or whether there is a smooth transi- 
tion between the two. There is a luminous tail to the 
Type Ic distribution: for example SN2010ay reached a 
peak luminosity of —20.2 mag, though still had a light 
curve consistent with being powered by nickel decay and 
did not show the spec troscopic features typical of SLSNe 
(|Sanders et al.l I2012D . If the dearth of intermediate- 
luminosity objects represents a real cutoff rather than 
a selection effect, this places constraints on any pro- 
posed mechanism for powering the SLSNe. Such a low- 
luminosity cutoff is not predicted by theoretical models 
of SLSNe; for example th e magnetar models presented in 
iKasen fc Bildster] ([2010D can reproduce a wide range of 
luminosities and timescales. 

We also note that PSl-lObzj only has a few, early 
light curve points that are brighter than —21 mag, but is 
clearly a spectroscopic member of the class of 2005ap-like 
hydrogen-poor SLSNe. This suggests that a definition 
based on a luminosity cut, as was suggested in iGal-Yaml 
( 2012), is artificial and that this class of objects is better 



— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — , 


1 1 1 1 1 1 1 1 1 1 ; 




























PS1-10bzj 























-21.0 



-21.5 -22.0 -22.5 

Peak absolute magnitude 



-23.0 



Figure 12. Distribution of peak observed absolute magnitudes 
for the hydrogen-poor SLSNe published to date. SN 2006oz is not 
included, as it was only observed on the rise and so only a lower 
limit < —21.5 mag is available. 

defined by spectroscopic features. 

5.2. The Host Galaxy Environments 

Of the 10 2005ap-like hydrogen-poor SLSNe published 
prior to this work, only five have detected host galax- 
ies and the uppe r limits on the undetected ones are 
M r > -18 mag (jNeill et al.l 120111) . It has been spec- 
ulated that this preference for low-luminosity environ- 
ments is really a preference (and perhaps requirement) 
for low-metallicity environments. This is supported by 
the host of SN 2010gx, the first SLSN host galaxy with 
a direct metallici ty measurement, with Z = 0.06Z Q 
(jChen et al.ll2013D h The low metallicity of 0.13 Z for 
the host of PSl-lObzj follows the same trend. 

To put these galaxy measurements in context, in 
Figure [T2] we plot different properties of the two 
SLSN host galaxies, compared to other galaxy sam- 
ples. The top left panel shows a mass-metallicity (M- 
Z) pl ot, with metallicity m e asure d by the R23 method 
from IKobulnickv fc Kewlevl (|2004D to facilitate compar- 
ison to different sample s, including core-collapse SNe 
ijKellv fc Kirshnerl 120121) and gamma - ray burst (GRB) 
host galaxies dLevesaue et al.1 12010bt iLeibler fc Bergerl 
I2010t iModiaz et al.ll2008l) . We plot here also the host 
of th e superluminous S N2007bi, with R23 metallicity 
from I Young et alJ (|2010l ) , and mass we estimated from 
the SDSS photometry of this host. It is worth noting 
that all three SLSN hosts have very similar R23 metallic- 
ities. They are all less massive and more metal-poor than 
the core-collapse SN hosts, but occupy a similar region 
as the gamma-ray burst (GRB) host galaxies. One im- 
portant caveat here is that the core-collapse SN hosts are 
generally at low redshift and contain a mix of hosts from 
targeted and untargeted surveys, so we would not nec- 
essarily expect them to follow the same mass-metallicity 
relation. 

The remaining three panels plot metallicity as mea- 
sured by the T e method, against either luminosity, 
mass, or the combination of mass and star formation 
rate that minimizes scatter in m etallicity (the so-called 
Fundamental Relation o r FMR; iMannucci et al.1 120101 : 
I Andrews fc MartinillMl ) . The host of PSl-lObzj is con- 
sistent with each of the nearby relations within its un- 
certainties, indicating that it is not an unusually metal- 
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poor galaxy given its mass, luminosity and star forma- 
tion rate. The similarity to GRB hosts may indicate that 
the two phenomena happen in similar environments, but 
the sample sizes here are small. The host of SN 2010gx 
stands out as more extreme than the host of PSl-lObzj 
in terms of metallicity, and falls below the nearby/SDSS 
relations in each case. The two SLSN hosts are the most 
separated on the FMR plot, due to the larger sSFR of 
the PSl-lObzj host. As such, the most striking common 
factor between the two galaxies is their low metallicities. 

We note that if low metallicity is an important fac- 
tor in producing this type of SLSN, this may present a 
challenge for models where the luminosity of the SN is 
powered by interaction with a dense wind. In particular, 
the mass loss would be unlikely to be driven by metal- 
line winds and so a different mass-loss mechanism would 
be required. 

6. CONCLUSIONS 

We show that PSl-lObzj is a hydrogen-poor superlumi- 
nous superno va, spectroscop i cally simila r to the objects 
descri bed in iQuimbv et al.1 (|2011| ) and iChomiuk et al.1 
(|2011ft . Compared to previous events, it has the fastest 
rise time and lowest peak luminosity. From our re- 
constructed bolometric lightcurve, we estimate the to- 
tal energy radiated over the time period observed to be 
~ 3.5 x 10 50 ergs, and the bolometric magnitude at peak 
to be about —21.4 mag. A magnetar model can fit the 
observed light curve, velocities and temperatures. Pro- 
posed interaction scenarios for SLSNe can also match 
the observed energetics but would require at least ~ 3 
Mq of hydrogen-poor circumstellar material. The lack of 
intermediate-width lines in the spectra, like with other 
SLSNe, also speaks agains this model. A normal, Ni- 
decay Ic model would require Mni = 7 M Q and the ejecta 
composition to be > 80 % Ni by mass, so although PSl- 
lObzj shows less extreme energetics than other hydrogen- 
poor SLSNe, radioactive decay is unlikely to be the pri- 
mary energy source. 

Like SN2010gx and PTF09cnd, PSl-lObzj developed 
a number of spectral features after peak. Our model fits 
these features with intermediate-mass elements Mg, Ca 
and Si, and Fe is also likely. We do not have the spectro- 
scopic coverage to determine whether these features at 
late times evolved into a more typical Type Ic SN spec- 
trum, as was seen in the other two objects. However, 
PSl-lObzj is interesting in the comparison to Type Ic 
SNe in the sense that it extends the distribution of SLSNe 
towards lower luminosities. Continuing to map out the 
low-luminosity tail of the SLSN population will be neces- 
sary to determine whether the two classes represent truly 
distinct phenomena, or whether there is a smooth con- 
tinuum between them. If the timescalcs of PSl-lObzj are 
typical for the lower-luminosity objects, this may present 
a challenge for finding such events as they will not stand 
out photometrically as much as higher-luminosity events 
and will require spectroscopic confirmation. 

The host galaxy of PSl-lObzj is detected both in 
our PS1 template images and in catalogs covering the 
ECDF-S. Combining this photometry with emission line 
measurements, we find that the host is a low luminos- 
ity (Mb — —18 mag; L ~ 0.05L*), low metallicity 
(Z = 0.13Z Q ), low stellar mass (M* w 2.4 x 10 7 M ) 
galaxy. It is forming stars at a rate of ~ 2—3 M yr _1 , re- 



sulting in a high sSFR (100 Gyr" 1 ). Archival HST imag- 
ing further reveal the host to be compact, with a physical 
diameter < 800 pc. While the metallicity is not as low as 
the host galaxy of the superluminous SN2010gx, the dis- 
covery of a second low metallicity host galaxy supports 
the hypothesis that metallicity may be important in the 
progenitor channel of SLSNe. Compared to the host of 
SN2010gx, the host of PSl-lObzj has a higher star for- 
mation rate, and is generally consistent with the mass- 
met all icity relation for starforming galaxies at lower red - 
shifts (|Mannucci et all 120101 : lAndrews fe Martinill2013f) . 
Further increasing the sample of SLSNe with well-studied 
host galaxies will be necessary to assess whether this 
metallicity trend holds, and shed light on the nature of 
these extreme explosions. 
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Figure 13. A comparison of the host of PSl-lObzj to other SLSNe and galaxy samples. The blue star and square show the hosts of 
PSl-lObzj and SN 2010gx, respectively. The top left panel is a mass-metallicity (M-Z) diagram, here plotting metallicity as calculated by 
the i?23 method in Kobulnicky & Kcwlcv (2004), to facilitate comparison to a broader sample. The b lack lines show the SDSS M-Z relation 
Tremonti ct al. (2004), where metallicity has been converted to the KK04 scale using the relations in Kcwlcy & Ellison (2008). The orange 
circles show hosts of core-collapse supernova (any type), with metallicities similarly converted from the Tremonti ct al. (2004) scale. Red 
triangles show GRB hosts; connected points indicate a dual solution for either the mass or the metallicit y. The other three plots show 
metallicity measured by the direct method, comparing to either the nearby luminosit y-metallicity relation (Guseva ct al. 2009; top right) 
the SDSS M-Z relation as measured by the direct method on spectra stacked by mass I Andrews fc Martinill20ia botto mTeft), or the SDSS 
"fund amental relation" , plotting metallicity against a combination of stellar mass a nd SFR that minim izes scatter {Xndrews & Martinil 
120131 ; b ottom right). Add itional samples shown are DEEP2 star-forming galaxies (Hoyos ct al. 2005), nearby Blue Compact Galaxies 
(BCGs; Kcwlcv ct al. 2007) and "Green Pea" galaxies (Amorin ct al. 2012). The host of PSl-lObzj is consistent with nearby relations in 
each case, and similar to the GRB host galaxies with direct metallicity measurements. 
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21.57 ± 0.11 


Pbl 


55586.1 


+ 13.6 


z 


21.59 ± 0.22 


Gemini-b/GMOb 


55586.3 


+ 13.7 


zpi 


21.55 ± 0.08 


Pbl 


55589.2 


+15.5 


zpi 


21.57 ± 0.08 


Pbl 


55595.2 


+19.1 




21.84 ± 0.09 


Pbl 


55627.5 


+38.7 




22.32 ± 0.12 


Magellan/IMACb 



a In rest-frame days, relative to maximum light on MJD 55563.65 



Table 2 

Log of bpectroscopic Observations 



UT Date 


Epoch a 


Instrument 


Wavelength Range 


Slit 


Grating 


Filter 


Exp. time 


Mean 


(YYYY-MM-DD.D) 


(days) 




(A) 


(") 






(s) 


Airmass 


2011-01-13.2 


6.7 


LDbb3 


3540-9450 


0.75 


VPH-all 


none 


3900 


1.3 


2011-01-25.1 


13.9 


GMOb-b 


3320-6140 


1.0 


B600 


none 


2400 


1.4 


2011-01-28.1 


15.7 


GMOb-b 


5890-10100 


1.0 


R400 


OG515 


3000 


1.1 


2011-04-02.0 


54.5 


GMOb-S 


5530-9830 


1.0 


R400 


OG515 


900 


1.9 


2011-04-03.0 


55.1 


GMOb-b 


5530-9830 


1.0 


R400 


OG515 


1800 


1.9 



a In rest-frame days relative to maximum light on UT 2011-01-02.7. 
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Table 3 

PSl-lObzj Host Galaxy Photometry 



UT Date Filter AB Magnitude Telescope/Instrument 







9pi 


Z 1. 


V7 4- 
O f it 


u. 


"1 Q 
lo 


rol 






'11 


24. 


00 ± 


0. 


.12 


PSl 






ip i 


23. 


76 ± 


0. 


.10 


PS1 






zpi 


22. 


73 ± 


0. 


05 


PSl 






ypi 




> 21. 


7 




PSl 


2011-11-29 


9' 


24. 


37 ± 


0. 


.08 


Gemini-S/GMOS 


2011- 


10-21 


r' 


23. 


86 ± 





18 


Magellan/LDSS3 


2011-09-20 


%' 


23. 


12 ± 





.07 


Gemini-S/GMOS 


2012-07-19 


z 1 


23. 


67 ± 


0. 


.15 


Magellan/IMACS 


2012-12-04 


J 




> 23. 


8 




Magellan/FourStar 


2011- 


12-07 


K 




> 22. 


7 




Magellan/FourStar 






F606W 


24. 


13 ± 


0. 


.05 


HST/ACS a 
HST/ACS a 






F850LP 


23. 


63 ± 


0. 


.06 






1/38 


24. 


89 ± 





08 


ESO MPG 2.2m/WFI b 






U 


24. 


86 ± 


0. 


.04 


ESO MPG 2.2m/WFI b 
ESO MPG 2.2m/WFI b 






B 


24. 


45 ± 


0. 


.02 






V 


24. 


44 ± 





.02 


ESO MPG 2.2m/WFI b 






/,> 


24. 


22 ± 





02 


ESO MPfl 2 2m /WFT b 






i 


Zo. 


Zo It 


n 
u. 


nK 
,uo 


jCjovJ IVlr^Vjx Z.zm/ Vvrl 






/ 

z 


2.3. 


39 ± 





13 


/~1 TIT /~\ A It, T " TTh 

GllO 4m/Mosaic-il u 






IA427 


21. 


25 ± 





.07 


bubaru/buprimc-Cam c 






IA445 


24. 


49 ± 


0. 


.07 


bubaru/buprimc-Cam 1 - 






IA464 


24. 


59 ± 





.14 


bubaru/buprime-Carm 






IA484 


24. 


45 ± 





03 


bubaru/buprime-Cam 






T A Knz 
IrttJUO 


21 


^±z in 


n 
u. 


.Uu 


Subaru / Suprime-Cam c 






IA527 


24. 


53 ± 


0. 


03 


Subaru / Suprime-Cam c 






IA550 


24. 


41 ± 





.05 


Subaru/Suprimc-Cam c 






IA574 


24. 


42 ± 


0. 


.06 


Subaru / Suprime-Cam c 






IA598 


24. 


21 ± 


0. 


.02 


Subaru/Suprime-Cam c 






IA624 


2.3. 


84 ± 


0. 


.02 


Subaru/Suprime-Cam c 






IA651 


24. 


08 ± 


0. 


.02 


Subaru / Suprime-Cam c 






IA679 


24. 


59 ± 


0. 


03 


Subaru / Suprime-Cam c 






IA709 


24. 


64 ± 


0. 


12 


Subaru/Suprime-Cam c 






IA738 


24. 


58 ± 


0. 


01 


Subaru / Suprime-Cam c 






IA767 


24. 


62 ± 


0. 


.10 


Subaru / Suprime-Cam c 






IA797 


2.3. 


78 ± 


0. 


06 


Subaru/Suprime-Cam c 






IA827 


22. 


25 ± 


0. 


03 


Subaru / Suprime-Cam c 






IA856 


24. 


68 ± 


0. 


.16 


Subaru / Suprime-Cam c 



a Data from GEMS survey catalog ( Rix ct al. 2004) 

Data from GaBoDs survey catalog jTaylor et al.|[2009D 
c Data from MUSYC survey catalog ( C ardamone ct al. 2010) 



Table 4 

Host Galaxy Emission Line Fluxes 



Line 


Jan 13 


Flux (10 16 erg s 1 cm 2 ) 
Jan 28 


Apr 3 


[O III]A5007 


6.43 ± 0.18 


6.36 it 0.07 


6.20 ± 0.08 


[O III]A4959 


1.89 ± 0.14 


2.04 ± 0.05 


2.32 ± 0.06 


H/3 


1.05 ± 0.16 


1.04 ± 0.05 


0.94 ± 0.08 


[O III]A4363 




0.14 ± 0.04 




H 7 


0.47 ± 0.07 


0.53 ± 0.04 


0.43 ± 0.05 


H<5 


0.20 ± 0.06 


0.22 ± 0.04 


0.24 ± 0.04 


He + [Nc III]A3968 


0.28 ± 0.08 


0.27 ± 0.05 


0.28 ± 0.06 


HC 




0.19 ± 0.05 




[Ne III]A3869 


0.48 ± 0.08 


0.37 ± 0.04 


0.58 ± 0.06 


[O II]A3727 


0.97 ± 0.09 


0.83 ± 0.06 


1.07 ± 0.08 



